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Abstract 



We prove that the quantum stress tensor for a massless scalar field 

\^ . in two dimensional non-self-similar Tolman-Bondi dust collapse and 

CN I Vaidya radiation collapse models diverges on the Cauchy horizon, if 

the latter exists. The two dimensional model is obtained by sup- 

^ . pressing angular coordinates in the corresponding four dimensional 

•/^ ! spherical model. 
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Q : 1 Introduction 
o 

! If a classical model of gravitational collapse results in the formation of a 
naked singularity, quantum effects can be expected to play a significant role 
^ ■ during the final stages of the collapse. One way to study these effects is the 
quantization of test matter fields in the background spacetime provided by 
the collapsing classical matter. The semiclassical approximation is expected 
to be valid up to Planck scales. In particular, one is interested in the behavior 
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of the stress tensor of the quantized field in the approach to the Cauchy 
horizon. The divergence of the vacuum expectation value of the quantized 
stress tensor on the Cauchy horizon signals an instability of the horizon, and 
suggests that back-reaction will prevent the naked singularity from forming. 

Two well-known examples of formation of naked singularities are the 
spherical collapse of inhomogeneous dust (the Tolman-Bondi model) and the 
spherical collapse of null dust (the Vaidya model). It is known for both these 
models that for certain initial data the collapse ends in a black hole and for 
other initial data it ends in a naked singularity 0, The stress-tensor 
of a quantized scalar field on these background spacetimes has been investi- 
gated, by specializing to the case of 2-d self-similar collapse. The restriction 
to 2-d is similar to the geometric optics approximation in 4-d - the latter 
amounts to keeping only the / = modes. The 2-d spacetime is obtained by 
suppressing the angular coordinates in the 4-d spherical model. For a two 
dimensional model, explicit expressions for the vacuum expectation value of 
the stress-tensor can be obtained from the trace anomaly, by imposing con- 
servation of the stress tensor. The assumption of self-similarity allows double 
null coordinates to be constructed explicitly. Using these coordinates it has 
been shown that the outgoing quantum flux diverges on the Cauchy horizon, 
for self-similar Vaidya collapse 0, as well as for self-similar Tolman-Bondi 
collapse 0]. 

A priori, it may be the case that the divergence in the 2-d model could be 
because of the assumption of self-similarity. In this paper, we prove that this 
assumption can be relaxed, and that the outgoing quantum flux will diverge 
on the Cauchy horizon, for all initial conditions for which a naked singularity 
forms in the 2-d Vaidya and Tolman-Bondi models. 

The outline of the proof is as follows. Consider the collapse of a classical 
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non-self-similar spherical dust cloud and choose the initial conditions to be 
such that the collapse results in a naked singularity. We now construct 
a new initial distribution by replacing a spherical region by a self-similar 
distribution. The new distribution is hence a self-similar spherical region 
surrounded by part of the original distribution. The free parameter of the 
self-similar distribution is fixed by requiring the first and second fundamental 
forms to match at the boundary between the self-similar region and the 
original distribution. It is then shown that if the evolution of the original 
distribution results in a naked (covered) singularity, the evolution of the 
modified distribution also results in a naked (covered) singularity. We show 
that, in general, the density of the cloud will change discontinuously at the 
boundary, but this change will be finite, and not infinite. 

We next consider the quantum stress tensor for a massless scalar field 
on the classical background dust spacetime. As has been shown earlier, 
the stress tensor diverges on the Cauchy horizon for a self-similar model. 
Consider now the modified distribution (self-similar region surrounded by 
non-self-similar region) mentioned above. We show that the nature of the 
divergence in the self-similar region is such that it implies a divergence in 
the outer non-self-similar region as well. Then, by considering a family of 
initial distributions, the size of the self-similar region is shrunk to zero - for 
each distribution in the family there is a divergence on the Cauchy horizon. 
The limiting distribution, in which the self-similar region disappears entirely, 
is the original distribution, and this also has a divergence on the Cauchy 
horizon. 

The plan of the paper is as follows. In Section 2 we construct the modified 
distribution which includes a self-similar spherical region in the interior. In 
Section 3 we obtain the quantum stress tensor for a test scalar field on this 
modified distribution and prove that it diverges on the Cauchy horizon for a 
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non-self-similar dust cloud. 



2 The Modified Distribution 



In this Section, we show how the modified initial distribution is constructed, 
for the cases of marginally bound and non-marginally bound dust collapse, 
as well as for the Vaidya model. 

2.1 Marginally Bound Dust Collapse 

The collapse of a spherical dust cloud is described by the Tolman-Bondi 
line-element, using comoving coordinates {t,r,6,(j)), 

ds" = -dt" + -^-—dr'' + R\t)dn\ (1) 
1 + /(r) 

where R{t,r) is the area radius at time t of the shell labeled r, and /(r) is 
a free function, satisfying / > — 1. The marginally bound solution is one 
for which /(r) = 0. The only non-zero component of the energy- momentum 
tensor is the energy density p(t,r), which satisfies the Einstein equation 

F'(r) 

where F{r) is another free function, and has the interpretation of being the 
mass to the interior of the shell r. The only other Einstein equation is 

R' = ^ + nr). (3) 

Let us consider the collapse of marginally bound dust cloud, starting at 
a time and having an initial density distribution po{r), for r < r^, where 
Tfe is the boundary of the star. Integrating @ gives the solution 



R'^%r) = ^^Fir)it,ir)~t) (4) 



for the evolution of the area radius of the shell r. to{r) is a function of 
integration, to be determined by choosing an initial scaling, R{ti,r), at the 
start of collapse. The area shrinks to zero at t = tQ{r), resulting in the 
formation of a curvature singularity. It is the singularity at r = 0, the 
central singularity, which is of interest to us, as this has been shown to be 
naked for some initial conditions. Specifically, it has been shown [|l| that if 
the initial density distribution po{R) has a Taylor expansion 

Po{R) = Po + PiR+Ip2R^ + Ip3R^ + ... (5) 
2 

then the singularity is at least locally naked if one of the following conditions 
is satisfied: (i) pi < 0, or (ii) pi = 0,p2 < 0, or (iii) pi = P2 = 0,p3 < and 
^ = V^p^/Apl^"^ is less than or equal to —25.9904. The singularity may or 
may not be globally naked. We are interested in showing that in either case 
the outgoing quantum flux diverges on the Cauchy horizon. 

The self-similar solution, i.e. one for which the spacetime of the collapsing 
cloud possesses a homothetic Killing vector field, is a special case of the 
marginally bound solution. If we choose the scaling in such a way that 
^o(^) = ^, then the mass function of the self-similar solution is of the form 
Fss{r) = Ar, with A a non-negative constant. All dimensionless quantities 
are functions of t/r. The central singularity forms at t = 0, and is known to 
be (globally) naked for A < Ac = 0.1809 0. It can also be shown that the 
initial density distribution for a self-similar cloud is of the form 

po (i?) = po + Pai?' + P6^' + P9^' + • • • (6) 

In order to construct the modified initial density distribution, we start 
with the original distribution po{r) and replace a central region, up to some 
r = Tc, by the self-similar solution (with rc < rb). For Vc < r < Vb, the 
distribution continues to be the original distribution po{r). There is only one 
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free parameter in the self-similar solution, namely A, and this is determined 
by requiring that the total mass in the self-similar region is equal to the 
total mass contained in the original distribution, up to r = Vc- That is. 
Arc = F{rc). 

We would now like to ensure that the determined value of A is such that in 
the limit Vc going to zero, the modified distribution admits a naked singularity 
if and only if the original distribution po{r) admits a naked singularity. This 
is essential because we are interested in examining properties of the Cauchy 
horizon in the original distribution. Hence it is natural to demand that the 
modified distribution possess a Cauchy horizon if and only if the original 
distribution does. The limiting value Aq, of A, is clearly dF{r)/dr\r=o, and 
should satisfy the condition A < Ac if and only if the original distribution 
Po{r) admits a naked central singularity. 

To check this, we go back to the solution and assume, for simplicity, 
the scaling to be such that tQ{r) = r, so that the solution can be written as 

R' = lF{r) {t - rf . (7) 

We substitute this in (0) and assuming the collapse to begin at t = tj„, find 
the initial density, pQ{r), near the center, to be 

2F{r) 



Po[r) 



3t 



1 + 



F'{r) t, 



(8) 



We compare this with the form (|^) for the initial density, after using (|^) with 
t = tin in (i)- This comparison allows us to deduce the form of F{r) near 
r = 0, from which the limiting value Aq = dF{r)/dr\r=o can be worked out. 
We find that for two of the naked cases, (i) pi < 0, and (ii) pi = 0, p2 < 0, we 
get Aq = 0, which is a special case of the naked self-similar model. For the case 
(iii) pi = p2 = 0,p3 < we get Aq^^ = — 8po''^/-\/3p3, which implies that the 
chosen self-similar distribution is naked if and only if the original distribution 
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of type (iii) is naked. For the covered case, which is pi = P2 = Ps = 0, we 
get that Ao is infinite, which is a covered self-similar distribution. Hence it is 
shown that the modified distribution admits a naked singularity if and only 
if the original distribution does. 

Next, we discuss the question of the matching of the self-similar region 
and the non-self-similar region, at the boundary Vc- By comparing the metrics 
of the two spacetimes at this boundary we find that the area radii in the self- 
similar metric and in the non-self-similar one will be equal if the two mass 
functions are equal. Since the masses have been chosen to be equal, that 
ensures the matching of the first fundamental form (i.e. the line element on 
the boundary). The second fundamental form (i.e. the extrinsic curvature) 
is defined to be 

i^^. = -^(^.;/3 + Cft«)^>^ (9) 

where ^„ is a unit normal to the boundary, and h'^ = 6^ — ^"^/3 is the pro- 
jection tensor. The only non-zero component of the extrinsic curvature is 
Kg = = —1/R, which matches for the two metrics, since the mass func- 
tions have been chosen to be equal. Hence we are assured that the first and 
second fundamental forms match on the boundary. 

We note that the derivative F'{r) will in general not match at the bound- 
ary, when calculated for the self-similar region, and for the non-self-similar 
region. As a result, the density function p{t,r), given by (^, will in general 
be discontinuous at the boundary. However, it is important for our purposes 
to note that this discontinuity will be finite, since F'{r) will be finite. 

2.2 Non- marginally Bound Dust Collapse 

The metric function /(r) which appears in (|l|) is now non-zero. Equation (|^) 
can be solved exactly, for given F{r) and f{r). Whereas F{r) is determined 
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from the initial density distribution, /(r) gets determined from the velocity 



Given an initial distribution of density and velocity which is non-marginal, 
one cannot straightaway match it to a self-similar region at a boundary Tc, 
because for the former distribution we have /(r) ^ 0, and for the latter we 
have /(r) = 0. Matching of the line-elements requires that both F(r) and 
/(r) match on the boundary between the two regions. 

Thus, given the original non-marginal distribution {po{r), f{r)) we con- 
struct the modified distribution as follows. Replace the original model by a 
self-similar one, as before, for r < Vc- For some < r^, introduce a non- 
marginal dust distribution {po{r), fs{r)) in the region Tc < r < r*, which 
has the property that poir) (and hence F{r)) is the same as in the original 
distribution; /^(rc) = and /s(r*) = /(r*). As before, the self-similar model 
is chosen such that Arc = F{rc). The introduction of this sandwiched re- 
gion ensures that the first and second fundamental forms are matched at the 
boundaries Tc and r*. The limiting value of A, as Tc goes to zero, is again 
given by Aq = dF{r)/dr\r=o, and is calculated by letting both Tc and r* go 
to zero, while always keeping r* > Tc- By carrying out an analysis similar 
to the marginal case, it can be shown that Aq is such that the introduced 
self-similar model is naked if and only if the original distribution is naked. 



The collapse of a null dust cloud is described by the Vaidya metric 



where the mass function m(t') depends on the advanced time coordinate v. 
The mass function is zero for v < 0, and constant for v > Vb- Thus the cloud is 



distribution R at the onset of collapse. 



2.3 The Vaidya Model 




(10) 
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bounded in the region < v < Vb, and its exterior is Schwarzschild spacetime. 
A curvature singularity forms when the inner boundary of the cloud hits 
R = 0. This singularity is known to be naked for dmlv) / dv\v=o < 1/16 
[0] and covered for dm{v) / dv\y=o > 1/16. The special case of a self-similar 
model is described by the linear mass function, m(t>) = fxv, < v < Vb- 

Given a non-self-similar model m{v), we construct the modified distri- 
bution by replacing the original model by a self-similar one in the region 
< V < Vc, with Vc < Vb- The free parameter /i in the self-similar model is 
fixed by demanding fiVc = m{vc). The limiting value of /i, as Vc goes to zero, is 
dm{v)/dv\t,=Q, and it is apparent that the introduced self-similar distribution 
is naked if and only if the original distribution m{v) is naked. 

Since the boundary Vc is a null hypersurface, the matching of the space- 
times at the boundary is done by comparing, not the line elements, but the 
affine parameter for outgoing or ingoing null geodesies. We now show that 
the matching of the affine parameters on the boundary Vc implies that the 
two mass functions should be equal. For this purpose, it is convenient to re- 
strict to the 2-d line element obtained from (p!0| ) by suppressing the angular 
coordinates, and to write it in double null coordinates u, v, i.e. 

ds"^ = C'^{u,v) dudv. (11) 

The function C{u,v) satisfies the differential equation 

A _ 2m{v) \ 2m{v)_ C^v _ 

[ R )^ R^C^ ^ ^ ' 

By integrating the geodesic equation for the metric (|lTD it is shown that the 

affine parameter along ingoing and outgoing null geodesies is of the form 

P = a J C^du + b , 1 = ^1 C^dv + d. (13) 

Matching of the affine parameter at the boundary Vc hence requires that 
C^du = C^du, and the use of this equality in the differential equation (0) 

9 



gives the result that the mass functions of the two regions (self-similar and 
non-self-similar) must match. Further, it may be shown, as in the dust case, 
that the extrinsic curvatures match at the boundary between the two regions. 

We have now completed the construction of the modified classical distri- 
bution, which will be used to prove the divergence of the outgoing quantum 
fiux on the Cauchy horizon. 

3 The Quantum Stress Tensor 

We restrict attention to the two dimensional spacetime (Tolman-Bondi or 
Vaidya) obtained by suppressing angular coordinates in the four-dimensional 
spherical spacetime. The expectation value < Om|T^i/|Oj„ > of the energy- 
momentum tensor of a quantized scalar field in the Minkowski vacuum 1 j„ > 
can be calculated from the trace anomaly, and Wald's axioms. The trace 
anomaly is equal to 7?./247r, where TZ is the Ricci scalar for the background 
spacetime 0. The two dimensional spacetime can be expressed in terms of 
global null coordinates u and v as 



These coordinates are chosen such that the initial quantum state of the scalar 
field , which is the standard Minkowski vacuum |0j„ > on X~, is the vacuum 
with respect to the normal modes of the scalar wave equation in u, v coordi- 
nates. The components of < T^^ > are given by 



ds^ = {u, v) dudv. 



(14) 




(15) 




(16) 



(17) 
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Consider a situation where two solutions for the background metric are 
matched across a hypersurface, hke a geodesic, for instance, which in fact will 
be the case for our models. Let Cf and C| be the metrics in the two regions. 
In either of the regions where the solutions are given, the expressions (|15|-0) 
are analytic. Now consider the hypersurface (boundary) where the matching 
has taken place. Assuming the line element on the boundary and the second 
fundamental form of the boundary being calculated in both of the regions 
to be identical respectively, one obtains Ci = C2 = C at the boundary. The 
difference in the quantum stress tensor components can be expressed at the 
boundary ( e.g. from equation 1). 

< Tim > - < T2UU >= -l/l2nC (1/Ci - l/Cs),^,^ (18) 

Since both the metrics are analytic in their respective regions (we only require 
their second order partial derivatives to be finite) and if one extends each 
one of them smoothly across the boundary, one finds that the expression 
above is certainly finite. So, the discontinuous change in < T^i^i > across the 
boundary will be finite. The same can be said of the rest of the components. 

We next introduce double null coordinates in the various regions of the 
modified distribution under consideration. Consider first the case of marginally 
bound dust. The various regions are the introduced self-similar region (hence- 
forth labeled 2) in the coordinate range < r < r^, the original distribution 
in the region (henceforth labeled 1) < r < r^, and the Schwarzschild region 
(henceforth labeled 0) r > ri,. We write the line-elements in these regions as 

ds'^ = {U2, V2)dU2dV2 (19) 

in region 2, as 

ds'^ = B^{ui,vi)duidvi (20) 
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in region 1, and finally, as 



ds'^ = D^{u,v)dudv (21) 

in region 0. 

The relationship amongst these coordinates can be given, following the 
procedure described in Birrel and Davies 0. Because of matching the first 
fundamental form at various boundaries, we have 

ui = ai{u), v = Pi{vi), (22) 

and 

U2 = a2{ui), Vi = P2{V2). (23) 

The center of the cloud being given by both 

u = v (24) 

and 

U2 = V2- 2Ro (25) 

in a manner similar to Birrell and Davies, we obtain a relation between ii 
and Ml 

u = Pi iP2 (a2 (Ml) + 2i?o)) (26) 
as well as a relation between it and u 

u = Pi {132 («2 («i {u)) + 2Ro)) . (27) 

Starting from the expressions (|15|-|1^) for the expectation value of the 
stress-tensor, and using the coordinate relationships given above, we obtain 
the following expressions for the T^^, in various regions. In region 2, 

< T^,u, >= -Fu, {A^) + Fu, (J^^ , (28) 
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dv 



1 / du dv' 



/12 42 /12 

,U2,V2 ,U2 ,V2 



A" 



(29) 
(30) 



The function Fxiii) is defined as 



• Hi: 



(31) 



In region 1, 



1 I du dv 



2 r d2 o2 o2 

jUlj-Ul ,Ul ,Vl 

~B2 
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(32) 
(33) 
(34) 



Similarly in region 0, 



< T^u >= F^ [d^) + a'F^, {(3') + («') 

< >= -F, [d^) , 
< T„„ >= -l/247r (in (d^ 



(35) 

(36) 
(37) 



where ' indicates first derivative with respect to the argument, and 



«()=a2(«i()), /3()=A(/?2()). (38) 

Let us first consider region 2; we have the self-similar metric there accord- 
ing to the procedure adopted. Consider a parameter value A for which the 
collapse results in a naked central singularity. From Barve et al. we know 
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that the component < T„2„2 > diverges on the entire Cauchy horizon in the 
region (the other components remaining finite). It is obvious that < Tuu > 
will also diverge there. As pointed out above, it is important now to note 
that the discontinuity in the quantum stress tensor, at the boundary rc, is 
finite at the boundary between regions 2 and 1. This is because the various 
expectation values are determined by the trace anomaly and derivatives (up 
to second derivative) of the metric. The metric derivatives are finite on the 
Cauchy horizon, for r > 0. Hence the quantum stress tensor component 
< ^Miui > at the boundary, in the limit of approach from region 1, will be 
divergent. 

We now examine the expression for < T^^uj > , i.e. equation (^2]). The 
first term in this expression is finite, because we demand the metric com- 
ponent to be at least C^. This term is finite not only at the event on the 
boundary but all over the null ray (Cauchy horizon) emanating from that 
event (in fact, all over the region 1). 

Now, the second statement above implies that the second term in equation 
( p2D diverges at the intersection of the boundary with the null ray (Cauchy 
horizon). This term is a function of only the retarded null coordinate ui. 
Since it diverges at one event (on the boundary), it diverges all along the 
outgoing null ray {ui = constant). The first term, as mentioned before, is 
finite there. Hence the tensor component < Ty^^^-^^ > diverges all along the 
Cauchy horizon in region 1 as well. 

Finally, we consider a family of modified distributions, each with a suc- 
cessively smaller value of the boundary coordinate Tc, so that in the limiting 
case, Tc tends to zero. For each family in the distribution, there is a diver- 
gence of the outgoing fiux, in region 1. In the limit that Vc tends to zero, we 
recover the original non-self-similar distribution po{r) which admits a naked 
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singularity, and by virtue of the construction given here, has a divergence of 
< Tuj^ui > on the Cauchy horizon. 

We can make a similar argument at the second boundary, between region 
1 and the Schwarzschild region 0. Thus, the divergence occurs all over the 
Cauchy horizon, to whichever extent it exists. In particular, if we have a 
globally naked singularity, there is a divergence at the intersection of the 
Cauchy horizon with X"*". 

Also, we note that this argument works for any number of regions with 
different metric solutions matched at the boundaries. This is important 
from the point of view of the extra shell we needed to introduce in the non- 
marginally bound Tolman Bondi case. The finiteness of the rest of the tensor 
components can be argued on similar lines. 

As regards the Vaidya metric, the entire argument is the same except for 
the fact that the boundaries between the regions are ingoing null geodesies. 
In fact, PO being the identity function the calculation is simplified con- 
siderably. 

4 Conclusion 

We have shown that the quantum stress tensor diverges on the Cauchy hori- 
zon, if it exists, in non-self-similar Tolman Bondi dust and Vaidya radiation 
collapse in two dimensions. There is no direct way of deducing this by analyt- 
ical calculations of the expressions in the general case. However, we employ 
a limiting process of approaching the required spacetime metric via patching 
up tractable solutions. From this technique, it appears that the divergence 
technically results from the metric rendered non-invertible at the Cauchy 
horizon in the self-similar portion. This does not seem to be the ultimate 
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reason for the divergence, for we see that the divergence persists even after 
the hmit to the actual metric is taken. Needless to say, the divergence does 
not seem to be ultimately a result of the self-similar nature although that 
does help in making the divergence evident in the calculations. 
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